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Abstract 

We report on a study of the isospin- violating and conserving decays of the J/ ip and tp' charmonium state 
to AA7T and AA77, respectively. The data are based on 225 million J/ip and 106 million if}' events that 
were collected with the BESIII detector. The most accurate measurement of the branching fraction of the 
isospin- violating process J/tp — > AA71- is obtained, and the isospin-conserving processes J/ip — >■ AA77 
and ip' — > AAt] are observed for the first time. The branching fractions are measured to be B(J/ip — > 
AAvr ) = (3.78 ± 0.27 stat ) ± 0.29 sys ) x 1(T 5 , £(J/Y> -> AA77) = (15.7 ± 0.79 stat ± 1.52 sys ) x lO" 5 
and B(ip' ->■ AA77) = (2.47 ± 0.34 stat ± 0.19 sys ) x 10~ 5 . No significant signal events are observed for 
ip' — > AAtt decay resulting in an upper limit of the branching fraction of B(ip' — > AAvr ) < 0.29 x 10~ 5 
at the 90% confidence level. The two-body decay of J/ip — > S(1385)°A + c.c. is searched for, and the 
upper limit is B(J/if> -»• S(1385)°A + c.c.) < 0.81 x 10~ 5 at the 90% confidence level. 

PACS numbers: 13.25.Gv, 12.38.Qk, 14.20.Gk 
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I. INTRODUCTION 



The charmonium vector meson, J /if), is usually interpreted as an SU(3) singlet cc bound 
states with an isospin 1=0. Systematic measurements of its decay rates into final states that 
are isospin violating are of particular interest, since these results will provide a sensitive probe 
to study symmetry-breaking effects in a controlled environment. In this paper, we present a 
systematic study of isospin-conserving and violating decays of charmonium vector mesons into 
baryonic decays accompanied by a light pseudoscalar meson, namely J/if)(if>') — > AA77 and 
J /if) (if)') — > AA7T , respectively. 

This work is for a large part motivated by a controversial observation that was made in the past 
while studying the baryonic decay of the J/if>. Surprisingly, the average branching fraction of the 
isospin violating decay of J /if) — > AAir measured by DM2 [ 1 ] and by BESI |2J] was determined to 
be B(J/if> — > AA71- ) = (2.2 ±0.6) x 1CT 4 , while the isospin conserving decay mode J /if) — > AA77 
was not reported by either experiment. In 2007, the decays of J /if) and if)' to the final states with 
a AA pair plus a neutral pseudoscalar meson were studied using 58 million J/ip and 14 million 
if)' events collected with the BESII detector [3]. The new measurement suggested that the two 
previous studies of J /if) — > AAti may have overlooked the sizable background contribution from 
J /if) — > S 7T°A + cc. The BESII experiment removed this type of background contribution 
and only a few statistically insignificant J /if) — >■ AA71- signal events remained, resulting in an 
upper limit of B(J/if> — > AAvr ) < 0.64 x 10~ 4 . Moreover, the isospin conserving decay mode, 
J /if) — >■ AAi], was observed for the first time with a significance of 4.8cx. However, signal events 
of the channels if)' — > AAn° and if)' — > AA77 were not observed by BESII, and resulted in upper 
limits of B(ip' AAvr ) < 4.9 x 10~ 5 and B(if)' AA77) < 1.2 x 10" 4 . 

In 2009, BESIII collected 225 million J/ip |4|] and 106 million if)' |5j] events. These samples 
provide a unique opportunity to revisit these isospin conserving and violating decays with im- 
proved sensitivity to confirm the previous observations in J /if) decays with BESII. The ambition 
is to investigate as well the same final states in if)' decays with the new record in statistics, and look 
for possible anomalies. A measurement of these branching fractions would be a test of the "12%" 
rule The data allow in addition a search for the two-body decays J /if) — > S(1385)°A + cc. 
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II. EXPERIMENTAL DETAILS 



BEPCII is a double-ring e + e~ collider that has reached a peak luminosity of about 0.6 x 
10 33 cm _2 s _1 at the center of mass energy of 3.77 GeV. The cylindrical core of the BESIII de- 
tector consists of a helium-based main drift chamber (MDC), a plastic scintillator time-of-flight 
system (TOF), and a CsI(Tl) electromagnetic calorimeter (EMC), which are all enclosed in a su- 
perconducting solenoidal magnet providing a 1.0 T magnetic field. The solenoid is supported 
by an octagonal flux-return yoke with resistive plate counter muon identifier modules interleaved 
with steel. The acceptance for charged particles and photons is 93% over An stereo angle, and the 
charged-particle momentum and photon energy resolutions at 1 GeV are 0.5% and 2.5%, respec- 
tively. The detector is described in more detail in [|70. 

The optimization of the event selection criteria and the estimates of physics background sources 
are performed through Monte Carlo (MC) simulations. The BESIII detector is modeled with the 
GEANT4 toolkit [|8|,|9j]. Signal events are generated according to a uniform phase-space distribution. 
Inclusive J/tp and ip' decays are simulated with the kkmc flp| generator. Known decays are 
modeled by the evtgen 111 ill generator according to the branching fractions provided by the 
Particle Data Group (PDG) 13, and the remaining unknown decay modes are generated with the 



lundcharm model [131. 



III. EVENT SELECTION 

The decay channels investigated in this paper are J ftp (ip') — » AA7r° and J/ip (tp') — > AAr/. 
The final states include A, A and one neutral pseudoscalar meson (7r° or rf), where A (A) decays to 
n~p (n + p), while the 7r° and i] decay to 77. Candidate events are required to satisfy the following 
common selection criteria: 

1 . Only events with at least two positively charged and two negatively charged tracks are kept. 
No requirements are made on the impact parameters of the charged tracks as the tracks are 
supposed to originate from secondary vertices. 

2. The transverse momenta of the proton and anti-proton are required to be larger than 
0.2 GeV/c. Tracks with smaller transverse momenta are removed since the MC simulation 
fails to describe such extremely soft tracks. 
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3. Photon candidates are identified from the reconstructed showers in the EMC. Photon ener- 
gies are required to be larger than 25 MeV in the EMC barrel region (| cos6>| < 0.8) and 
larger than 50 MeV in the EMC end-cap (0.86 < | cos#| < 0.92). The overlapping showers 
between the barrel and end-cap (0.8 < | cos6>| < 0.86) are poorly reconstructed, there- 
fore, excluded from the analysis. In addition, timing requirements are imposed on photon 
candidates to suppress electronic noise and energy deposits from uncorrelated events. 

4. The A and A candidates are identified by a reconstruction of decay vertices from pairs of 



oppositely charged tracks pn~ and pn + IU4H . At least one pn~ and one pn + candidate are 
required to pass the A (A) vertex fit successfully by looping over all the combinations of 
positive and negative charged tracks. In the case of multiple A A pair candidates, the one 
with the minimum value of (M p7T - — M A ) 2 + (M pn + — M A ) 2 is chosen, where .\ /\(.U A ) is 
the nominal mass of A(A), obtained from the PDG II 1211 . 

5. To further reduce the background and to improve the resolution of the reconstructed parti- 
cle momenta, candidate signal events are subjected to a four constraint energy-momentum 
conservation (4C) kinematic fit under the hypothesis of J/ip (ip') — > AA77. In the case of 
several combinations due to additional photons, the one with the best xlc vame is chosen. 
In addition, a selection is made on the xlc- ^ ts vame is determined by optimizing the signal 



significance S/yS + B, where 5 (B) is the number of signal (background) events in the 
signal region. This requirement is effective against background with one or several addi- 
tional photons like J/ip, ip' -» S°7r°A + c.c. (S° ->■ 7A) or J/ip, ip' ->■ AA + (n > 4) 
decays (for instance J/ip, ip 1 £(1385)°£(1385)°, etc.). For J/if) AAtt , 

backgrounds are suppressed by requiring xlc < 40 (see Fig. HJa)). For J /if) — > AAr/, 
the requirement is set to xlc < 70 (see Fig. djb)). For if>' — > AA71- , due to the peaking 
background if>' — > S°7r°A + c.c. the xlc ^ s required to be less than 15 (see Fig. [Etc)). For 
if)' — > AAr], we select events with xlc < ^ ( see Fig-Ed)). 

Followed by the common selection criteria, a further background reduction is obtained by ap- 
plying various mass constraints depending on the channel of interest. To select a clean sample of 
A and A signal events, the invariant masses of pn~ and pn + are required to be within the mass 
window of \M piT — M\\ < 5 MeV/c 2 . Here, the invariant mass is reconstructed with improved 
momenta from the 4C kinematic fit. The mass resolutions of A and A are about 1.0 MeV/c 2 . For 
J/ip — > AA71- , a mass selection of |M p7r o — 1 189.0 1 > 10 MeV/c 2 is used to exclude background 
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from J /if} — > £ + 7r~A + c.c. (S + — > pn°) which can form a peak near the n° mass. The back- 
ground from J /if) is removed by selecting events with M AA < 2.8 GeV/c 2 as shown in 
Fig.[2£a). For J /if) — > AAr/, a selection of events with M AA < 2.6 GeV/c 2 rejects all background 
contributions from J/ip — > decays as shown in Fig.^b). For if)' — > AA7r° and ^' — >■ AAr/, 
events must satisfy the condition \M^° l l — 3097 1 > 8 MeV/c 2 to remove the background from 
if} 1 — > n + n~J/ip (J /ip — > ppn° and ppr/). The background from if)' — >■ jjj/if} {J /if) — ?• AA) and 
■0' — >■ £°£° is rejected by the requirement M AA < 3.08 GeV/c 2 . The AA invariant-mass distribu- 
tions for data and MC events from if)' — > AAvr , AA77, and E°E° are shown in Fig. [3] The scatter 
plot of M pw - versus M pn + after applying all selection criteria is shown in Fig. HI No visible signal 
of if}' — > AAtt is observed. 




FIG. 1 : The xlc distributions of 4C fits. Dots with error bars denote data, and the histograms correspond to 
the result of MC simulations, (a) J /if) — > AA7r°. The dashed line is the dominant background distribution 
from J /if) — > E°7r°A + c.c. with MC simulated events, the arrow denotes the selection of x\c <40. (b) 
J /if) — > AAr/, the arrow denotes the selection of x\c <70. (c) if)' — > AA7r°. The dashed line is the 
dominant background distribution from if)' — > S°7r°A + c.c. with MC simulated events, the arrow denotes 
the selection of xlc (d) if)' — > AAr/, and the arrow denotes the selection of xlc <40. 
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FIG. 2: The AA invariant-mass, M AA , distributions for J/ip AA77 candidates. Dots with errors denote 
data. The dashed-line shows the result of MC simulated events of J /if) — > which is normalized 

according to the branching fraction from the PDG. (a) Histogram shows the MC simulated events of J ftp — > 
AA7r , where the arrow denotes the selection of M AA <2.8 GeV/c 2 . (b) Histogram shows the MC simulated 
events of J/ip — > AA77, and the arrow shows the selection of M AA <2.6 GeV/c 2 . 



IV. BACKGROUND STUDY 



Backgrounds that have the same final states as the signal channels such as J /if), if)' — > 
£ + 7r~A + c.c. are either suppressed to a negligible level or completely removed. Back- 
ground channels that contain one or more photons than the signal channels like J /if), if>' — > 
E(1385)°S(1385)°, have very few events passing event selection. The line shape of the 
peaking background sources, J /if), if)' — > £°7r°A + c.c, is used in the fitting procedure to estimate 
their contributions. The contribution of remaining backgrounds from non-AA decays including 
J /if), if)' — > n + n~ppn° (77) is estimated using sideband studies as illustrated in Fig. [4j The square 
with a width of 10 MeV/c 2 around the nominal mass of the A and A is taken as the signal region. 
The eight squares surrounding the signal region are taken as sideband regions. The area of all the 
squares is equal. The sum of events in the sideband squares, Wideband region* times a normaliza- 
tion factor / is taken as the background contribution in the signal region. The normalization factor 
/ is defined as 

j- -^signal region 

-^sideband region 

The normalization factor is obtained from phase-space MC simulations of J /if) [if)') — > 
pp% + n~7r° or ppir + n~ri with Af signa i reg ion as the number of MC events in the signal region and 
^sideband region as the sum of MC events in the sideband regions. 

With 44 pb 1 of data collected at a center-of-mass energy of E cm = 3.65 GeV, the contribution 
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FIG. 3: The AA invariant-mass, M AA , distributions for ip' — > AA77 candidates, (a) MC simulated events 
of ip' ->■ AA7T , (b) MC simulated events of ip' -> A A??, (c) MC simulated events of ip' -> and (d) 

data. The arrow denotes the selection of M AA <3.08 GeV/c 2 . The peak around the J/ip mass is from the 
decay of ip' — > jjj/ip, J/tp — > AA. 

from the continuum background is determined. From this data sample, no events survive in the 7r° 
or T) mass region in the two-photon invariant-mass, M 77 , distribution after applying all selection 
criteria. Therefore, we neglect this background. 

V. SIGNAL YIELDS AND DALITZ ANALYSES 

The 77 invariant-mass spectra of J/ip — > AA71- , AAi], ip' — > AA%° and AAr/ of the remaining 
events after the previously described signal selection procedure are shown in Fig. [5J A clear n° 
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FIG. 4: A scatter plot of M p7T + versus M pn - for J/ip and ip' data, (a) J/i/> ->■ AA7r°, (b) -> AAr], (c) 
V?' -)■ AA7T , and (d) ?// -)■ AAt/. 

and r\ signal can be observed in the J/ij) data. The if)' data set shows a significant r] signal, but 
lacks a pronounced peak near the n° mass. 

The number of signal events are extracted by fitting the M 77 distributions with the parame- 
terized signal shape from MC simulations. For J/tp (ip') — > AA7r°, the dominant peaking back- 
grounds from J/ip (ip') — > £°7r°A + c.c. are estimated by MC simulation. The fit also accounts 
for background estimates from a normalized sideband analysis. Other background sources are de- 
scribed by a Chebychev polynomial for all channels except ip' — > AAn° where there are too few 
events surviving. The fit yields 323 ± 23 7r° events, 454 ±23 t] events in J/ip data and 60.4 ± 8.4 r] 
events in ip' data. For ip' — > AA7r°, the upper limit on N w o is 9 at the 90% confidence level (C.L.) 
and is determined with a Bayesian method tl5|]. For ip' — > AArj, the change in log likelihood value 
in the fit with and without the signal function is used to determine the r] signal significance, which 
is estimated to be 10.5cx. 

To study the existence of intermediate resonance states in the decay of J/ 'ip — > AA7T , 
J/ip — »■ AAt] and ip' — > AAi] and to validate the phase-space assumption that was used in the MC 
simulations, we have performed a Dalitz plot analysis of the invariant masses involved in the three- 
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FIG. 5: The two-photon invariant-mass, M 77 , distributions in the 7r° and g mass regions for the channels (a) 
J/tp -> AA7r°, (b) J/-0 -> AA??, (c) -0' -> AA7T , and (d) ip' -> AAr?. Dots with error bars are data. The 
solid lines are the fit to data, and the dot-dashed lines are the signal shape determined from MC simulations. 
The hatched histograms are the background contributions obtained from a normalized sideband analysis. 
The dashed lines in J/ip, ip ~~ > AAvr correspond to the peaking background from S°7r°A. The long 
dashed lines denote other background contributions which are described by Chebychev polynomials. 



body decay. These results are shown in Fig. [6] For these plots, n° and t] candidates are selected 
within mass windows of 0.12 GeV/c 2 < M 77 <0.14 GeV/c 2 and 0.532 GeV/c 2 < M 77 < 0.562 
GeV/c 2 , respectively. In all the Dalitz plots, no clear structures are observed. A x 2 test is per- 
formed to confirm the consistency between data and the phase-space distributed MC events. The 
X 2 is determined as follows: 



X 



( n data _ n MC/ g f 
/ j „ data 



n MC \ ^data/MC 



, n, ' refers to the number 



where g is the scaling factor between data and MC (g 
of data/MC events in a particular bin in the Dalitz plot, and the sum runs over all bins. We 
divide the Dalitz plots into 8 bins. Boxes with very few events are combined into an adjacent bin. 
The x 2 /n.d.f. are equal to 1.1 and 2.1 for J/ip — > AAtt and J/ip — > AArj, respectively, which 
validates the usage of a phase-space assumption in the MC simulations. 



12 



> * 

CD 

CD 3.5 


f (a) 








> 

CD 

CD 


4 
3.5 


r (b) I «^gL_. I 












3 












3 


- ^JBBBBfe-* — 


<^< 2.5 












2.5 


2 












2 


= — tWHSPt — 


1.5 












1.5 


~1 , , , i l. < < , i , < < i , , < < , , , , l , 



2.5 3 



(GeV 2 ) 



4° 



(GeV 2 



> 

CD 

CD 



3 
2.8 











, — »* 

CvJ 

^ 38 

CD 3.6 
3.4 
er 3.2 

■ST 3 
2.8 
2.6 


















- ■» 

: ->m 


















111 
















■ 







2.6 2.8 



3.2 3.4 3.6 



3.8 4 

(GeV 2 ) 



2.6 2.8 



3.2 3.4 3.6 



3.8 4 

/2\ 



(GeV^ 



> 

CD 
CD 



r (e) 
























=, , , , 1 , , , , 1 


rt.lt.,,!.., 


, T T , , T T 1 , T 1 





2.5 



3 3.5 4 4.5 5 5.5 6 6.5 7 
«2 /n„\/2\ 



M' 



'All 



(GeV^) 



> 


6.5 


CD 




CD 


6 




5.5 




5 




4.5 




4 




3.5 




3 




2.5 



t (f),« 

= ^ 




... 




1 :« 




; ■/ - 

<M 




III 
■ 1 


jyilif 
SIP 


fig 






WMffi 


«f§i 





2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 

/2\ 



(GeV^) 



FIG. 6: Dalitz plots of the invariant masses M^ , , versus M^ ^ for the channels (a) J/ip — > AAtt 
(data), (b) J/i/> -> AAvr (MC), (c) J/> -»■ AAr? (data), (d) J/^ -> AAr? (MC), (e) ^' -> AAr? (data), and 
(f) V' -> AAr? (MC). See text for more details. 

We have studied the branching fraction of the decay J/tp — > S(1385)°A + c.c. by com- 
bining and analyzing the invariant-mass spectra of Air and An pairs as depicted in Fig. [7J 
For this analysis, ir° events are selected by applying a two-photon invariant-mass selection of 
0.12 GeV/c 2 < M 77 <0.14 GeV/c 2 . For the fit, the signal function is taken from a MC simula- 
tion of J/if> — > S(1385)°A + c.c, and the background function is taken from a MC simulation 
of J /if} — > AA7T . A Bayesian analysis gives an upper limit on the number of S(1385)°A + c.c. 
events of 37 at the 90% C.L.. 
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FIG. 7: A search for S(1385)° events by fitting the combined M^ w o and M A7r o invariant-mass distributions. 
Dots with error bars are data. The solid line is the fit to data. The hatched histogram is the signal func- 
tion obtained from a MC simulation, and the dashed line is the background function obtained from a MC 
simulation of J/ip — > AA7r°. 

VI. SYSTEMATIC ERRORS 

To estimate the systematic errors in the measured branching fractions of the channels of interest, 
we include uncertainties in the efficiency determination of charged and photon tracks, in the vertex 
and 4C kinematic fits, in the selection criteria for the signal and sideband region, and in the fit 
range. The uncertainties in the total number of J ftp and ip' events and in the branching fractions 
of intermediate state are considered as well. Below we discuss briefly the analysis that is used to 
determine the various sources of systematic uncertainties. 

• Tracking efficiency. We estimate this type of systematic uncertainty by taking the difference 
between the tracking efficiency obtained via a control channel from data with the efficiency 
obtained from MC simulations. The control sample J/ip — > pK~A + c.c. is employed to 
study the systematic error of the tracking efficiency from the A (A) decay. For example, 
to determine the tracking efficiency of the ir + tracks, we select events with at least three 
charged tracks, the proton, kaon and anti-proton. The total number of ir + tracks, iV°+, 
can be determined by fitting the recoiling mass distribution of the pK~p system, M^ coi f. 
In addition, one obtains the number of detected n + tracks, iV*+, by fitting M^ coi f, after 
requiring all four charged tracks be reconstructed. The n + tracking efficiency is simply 

N 1 

evr+ = j^r~- Similarly, we obtained the tracking efficiencies for n~ , p, and p. With 225 x 
10 6 inclusive MC events, we obtained the corresponding tracking efficiency for the MC 
simulation. The tracking efficiency difference between data and MC simulation is about 
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1.0% for each pion track. This difference is also about 1.0% for a proton (anti-proton) if its 
transverse momentum, P t , is larger than 0.3 GeV/c. The difference increases to about 10% 
for the range 0.2 GeV/c<P t <0.3 GeV/c. Conservatively, we take a systematic error due 
to tracking of 1% for each pion. For the proton (anti-proton), we use weighted systematic 
errors, namely 2% in J /if) -> AA7r°, 3.5% in J /if) ->■ AAr], 1.5% in if) 1 ->■ AA7T , and 2% in 
if)' ->■ AAr?. 

• Vertex fit. The uncertainties due to the A and A vertex fits are determined to be 1.0% for 
each by using the same control samples and a similar procedure as described for the tracking 
efficiency. 

• Photon efficiency. The photon detection efficiency was studied by comparing the photon 
efficiency between MC simulation and the control sample J /if) — > p°n°. The relative effi- 
ciency difference is about 1% for each photon Rial , which value was used as a systematic 
uncertainty. 

• Efficiency of the kinematic fit. The control sample of J /if) -> S° (E°) ->■ 7A (A) 
is used to study the efficiency of the 4C kinematic fit since its final state is the same as our 
signal. The event selection criteria for charged tracks and photons and the reconstruction 
of A (A) are the same as in our analysis. If there are more than two photon candidates 
in an event, we loop over all possible combinations and keep the one with the smallest 
value for (M 7 a — M 2 o) 2 + (M 7 ^ — Mgo) 2 . Furthermore, the remaining backgrounds are 
suppressed by limiting the momentum windows of S° and S°, i.e., \P^o — 980 1 < 40 MeV/c 
and |Ps° — 980 1 < 40 MeV/c. Figure [8] shows the scatter plot of M 7A versus for the 
inclusive MC events and J/ip data after applying all event selection criteria. The square in 
the center with a width of 10 MeV/c 2 is taken as the signal region. Almost no background 
is found according to the topology analysis from inclusive MC events. The candidate signal 
events for both data and MC events are subjected to the same 4C kinematic fit as that in 
our analysis. The efficiency of the 4C kinematic fit is defined as the ratio of the number of 
signal events with and without a 4C kinematic fit. A correction factor, he can be obtained 
by comparing the efficiency of the 4C kinematic fit between data and MC simulation, i.e., 
he — "ifer- The efficiency corrections corresponding to \ < 15, 40, 70 are (90.3 ± 0.8)%, 
(97.5 ± 0.6)% and (98.7 ± 0.3)%, respectively. The errors in the efficiency corrections are 
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taken as a systematic uncertainty. 
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FIG. 8: The scatter plot of M 7 a versus for J/tp -)■ S°E° for (left) data and (right) inclusive MC 
events. 



• Fit range. The it , 77, and S(1385)° yields are obtained by fitting the data around the corre- 
sponding mass value. By changing the mass ranges for the fits, the number of signal events 
changes slightly. These differences are taken as the errors due to the uncertainty of the fit 
range. 

• Signal and sideband regions. By changing the signal and sideband region from 5 MeV/c 2 x 5 
MeV/c 2 to 6 MeV/c 2 x6 MeV/c 2 , the number of fitted n°,r] and S(1385)° events changes 
slightly for data and MC. The differences in yield between the two region sizes are taken as 
systematic errors. 

• Background shape. A part of the background depicted in Fig. [5] is estimated by a fit with 
a third-order Chebychev polynomial. The differences in signal yield with a background 
function that is changed to a second-order polynomial, are taken as a systematic error due 
to the uncertainty in the description of the background shape. 

• Total number of J/ip and ip' events. The total numbers of J/ijj and ip' events are obtained 
from inclusive hadronic J/ip and ip 1 decays with uncertainties of 1.24% |4J] and 0.81% |5], 
respectively. 

All the sources of systematic errors are summarized in Table U The total systematic error is calcu- 
lated as the quadratic sum of all individual terms. 
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TABLE I: Systematic errors in the measurements of the branching fractions (%). 



Source 


J/tp AAtt 


J/ip -> AAt? J/V> - 


-> S(1385)°A+c.c. 


i>' -> AAtt 


ip' -> AAt? 


Photon efficiency 


2.0 


2.0 


2.0 


2.0 


2.0 


Tracking efficiency 


6.0 


9.0 


4.0 


5.0 


6.0 


Vertex fit 


2.0 


2.0 


2.0 


2.0 


2.0 


Correction factor of 4C fit 


0.6 


0.3 


0.6 


0.8 


0.6 


Background function 


0.6 


0.2 


1.5 


negligible 


2.5 


Signal and sidebands 


3.6 


1.7 


negligible 


9.1 


2.0 


Fit range 


0.6 


0.4 


negligible 


negligible 


1.5 


B(A^ Trp) 


0.8 


0.8 


0.8 


0.8 


0.8 


B(P -> 77 ) 


negligible 


0.6 


negligible 


negligible 


0.6 


£(£(1385)° -> Avr°) 






1.7 






Nj/^ 


1.24 


1.24 


1.24 














0.81 


0.81 


Total 


7.8 


9.7 


5.6 


10.9 


7.6 



VII. RESULTS 



The branching fraction of J/ip (ip r ) — > X is determined by the relation 

N obs [j/ip <y) -> x -> y] 



£(J/V> (V) -»> *) 



x w wo • b (x -> y) • e [j/^ (^') -* x y] • / 4C ' 

and if the signal is not significant, the corresponding upper limit of the branching fraction is 
obtained by 

1 1 V W ) ) Njh {V) -B(X^Y)- e[J/V m ^X^Y}-f 4C - (1.0 - a sys .) ' 

where, X ob,s is the number of observed signal events or its upper limit Nffll, Y is the final state, 
X is the intermediate state, e is the detection efficiency, and a sys , is the systematic error. The 
branching fraction of X — > Y is taken from the PDG Ill2n . Table [H] lists the various numbers that 
were used in the calculation of the branching fractions. With these, we obtain 
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TABLE II: Numbers used in the calculations of the branching fractions. 



Channel J/np -> AAtt J/ip -> AA?y J/^ -> S(1385)°A+c.c. V' -> AAtt V' -> AA77 



Number of events N obs /N$£ 



Efficiency e (%) 

/4C (%) 

^Vj/V- (xlO 6 ) 
Jfy (xlO 6 ) 
B(A -> Trp) (%) 
£(£(1385)° -> Att°) (%) 
S(7r°,r,^ 77 ) (%) 



323 
9.65 
97.5 
225.3 

63.9 

98.8 



454 
8.10 
98.7 
225.3 

63.9 

39.4 



< 37 
6.22 
97.5 
225.3 

63.9 
87.5 



< 9 
8.95 
90.3 



60.4 
14.64 
97.5 



106.41 106.41 
63.9 63.9 

98.8 39.4 



B(J/ip -> AAtt ) = (3.78 ± 0.27 (stat.) ± 0.29 (sys.)) X 10" 5 , 

B(iff' -> AAtt ) < 0.29 x 10~ 5 , 

B(J/ip ->■ AAr?) = (15.7 ± 0.79 (stat.) ± 1.52 (sys.)) x 10~ 5 , 

B(ip' AAr/) = (2.47 ± 0.34 (stat.) ± 0.19 (sys.)) x 10~ 5 , 

B{J/ip -»■ S(1385)°A + c.c) < 0.81 x 10~ 5 . 

Here the upper limits correspond to the 90% C.L.. 

With these results, one can test whether the branching ratio between the tp' and J/tp decays to 



the same hadronic final state, Qh, is compatible with the expected 12% rule [6]. We find a Qh for 
the channels AA7r° and AA77 of 

B(ib' -> AAtt ) 

Qh = -wryri aXoT < 10 - 0% 

B(J/4> -> AA7T ) 



at the 90% C.L., and, 



The errors reflect a quadratic sum of the systematic and statistical error, whereby some of the 
common sources of systematic errors have been canceled. Clearly, the isospin-violated decay 
AA7r° is suppressed in ip' decays, while Qh for the isospin-allowed decay, AA77, agrees with the 
"12%" rule within about la. 
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VIII. SUMMARY 



This paper presents measurements of the branching fractions of the isospin-violating and 
isospin-conserving decays of the J/ip and ip' into AA7T and AA^, respectively. The results to- 
gether with the measurements from previous experiments are summarized in Table [TTTJ We note 
that the earlier measurements of the branching fraction of the decay AA7r° by BESI and DM2 
likely overlooked a sizeable background contribution in their analysis as supported by the BESII 
and BESIII results. Hence, we claim that we have observed for the first time the two processes, 
J /ip — y AAtt and ip' — > AAr/. Moreover, the branching fractions of the J/ip — > AAr] decay is 
measured with a drastically improved precision. Its central value is lower than the BESII mea- 
surement by about \.5a. The branching ratios of J j 'ip — > AA7r° and ip' — > AAr] are consistent 
with previous upper limits, and the upper limit of ip' — > AAir is significantly more stringent than 
the BESII measurement. The isospin-violating decay modes, J/ip — > AAir and ip' — > AA7T , are 
suppressed relative to the corresponding isospin-conserving decay modes into AArj, albeit only by 
a factor of 4 in the case of the J/ip decay. In addition, we search for the isospin violating decays 
of J /ip S(1385)°A + c.c. and no significant signal is observed. 

TABLE III: A comparison of the branching fractions of this work with the results of previous experiments (x 10 -5 ). 
The first error is statistical and the second one indicates the systematical uncertainty. 

Experiments B(J/tp ->■ AAtt ) B(J/ip -» AAr]) B{ip' -> AAtt ) B(tp' AAr;) 

This experiment 3.78 ± 0.27 ± 0.29 15.7 ± 0.79 ± 1.52 < 0.29 2.47 ± 0.34 ± 0.19 

BESII [3J < 6.4 26.2 ±6.0 ±4.4 < 4.9 < 12 

BESI[2J 23.0 ±7.0 ±8.0 

DM2 [IJ 22.0 ±5.0 ±5.0 
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